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ABSTRACT 

The  foundation  of  the  theory  of  functionally  graded  plates  with  simply  supported  edges,  under  a 
Friedlander  explosive  air-blast  are  developed  within  the  classical  plate  theory  (CPT).  Within  the  development 
of  the  theory,  the  two  constituent  phases,  ceramic  and  metal,  vary  across  the  wall  thickness  according  to  a 
prescribed  power  law.  The  theory  includes  the  geometrical  nonlinearities,  the  dynamic  effects,  compressive 
tensile  edge  loadings,  the  damping  effects,  and  thermal  effects.  The  static  and  dynamic  solutions  are  developed 
leveraging  the  use  of  a  stress  potential  with  the  extended-galerkin  method  and  the  runge-kutta  method. 
Validations  with  simpler  cases  within  the  specialized  literature  are  shown.  The  analaysis  focuses  on  how  to 
alleviate  the  unwanted  effects  of  large  deformations  through  proper  material  selection  and  the  proper  gradation 
of  the  constituent  phases  or  materials. 


1.  introduction 

During  combat  situations,  the  structure  of  army  military 
vehicles  may  have  to  structurally  endure  the  effects  of  blast 
loading.  Advances  in  FGM’s  which  combine  the  properties 
of  two  dissimilar  materials  has  been  a  motivating  factor  in 
viewing  these  types  of  materials  as  a  viable  alternative  to  the 
current  isotropic  metallic  structures  being  utilized  in  the  hull 
and  armor  plating.  FGM’s  are  microscopically 
inhomogeneous  with  thermo-mechanical  properties  which 
vary  smoothly  and  continuously  from  one  surface  to  another. 
These  graded  structures  allow  the  integration  of  dissimilar 
materials  like  ceramic  and  metals  that  combine  different  or 
even  incompatible  properties  such  as  hardness  and 
toughness. 

In  this  paper,  the  foundation  of  the  nonlinear  theory  of 
functionally  graded  plate-type  structures  under  an  explosive 
air-blast  is  developed.  An  approximate  solution 
methodology  for  the  intricate  nonlinear  boundary  value 
problem  is  devised,  and  results  that  are  likely  to  contribute  to 
a  better  understanding  of  the  structural  behavior  under  an 
explosive  blast  with  beneficial  implications  towards  their 
improved  design  and  exploitation  are  presented. 

It  is  the  author’s  intent  within  this  paper  to  fill  in  some 
major  gaps  currently  existing  within  the  specialized 
literature. 


2.  BASIC  ASSUMPTIONS  AND  PRELIMINARIES 

The  plate  mid  surface  is  referred  to  a  cartesian  orthogonal 
system  of  coordinates  (x,y,z),  while  z  is  the  thickness 
coordinate  measured  positive  in  the  upwards  direction  from 
the  mid-surface  of  the  plate  with  h  being  the  uniform  plate 
thickness  and  y  is  directed  perpendicular  to  the  x-axis  in  the 
plane  of  the  plate.  See  figure  1  below. 


Figure  1:  A  simply  supported  functionally  graded  plate 
Shown  in  2-D  under  an  explosive  blast. 


The  nonlinear  elastic  theory  of  FG  Plates  is  developed 
using  the  classical  plate  deformation  Theory  [6].  It  is  also 
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assumed  that  the  EG  plate  is  made-up  of  ceramic  and  metal 
phases  whose  material  properties  vary  smoothly  and 
continuously  across  the  wall  thickness.  By  applying  the  rule 
of  mixtures,  the  material  properties  such  as  Young’s 
Modulus,  Density,  and  Poisson’s  Ratio  are  assumed  to  vary 
across  the  wall  thickness  as 

P{z)=Pcyc{z)+Pnym{z),  (la) 

In  which  and  denote  the  temperature-dependent 
material  properties  of  the  ceramic  and  metallic  phases,  of  the 
plate,  respectively  and  may  be  expressed  as  a  function  of 
temperature  [7,  9]  as 

P=Pq(Pj7’‘‘  -i-1-i-Pi7’-i-P27’^+^37’^)  (lb) 

Where,  Pq,P_i,Pi,P2,  and  P3  are  the  coefficients  of 
temperature  T(K)  and  are  unique  to  the  constituent  materials. 

(z)  and  (z)  are  correspondingly,  the  volume  fractions 
of  the  ceramic  and  metal  respectively,  fulfilling  the  relation 

yc{z)+y,n{z)=i-  (3) 

By  virtue  of  (3),  Eq.  (la)  can  be  expressed  as 

P{z,  t)  =  [P^  {T)  -  P^  {T)W^  (z)  +  P^  {T) .  (4) 

By  observation,  one  can  deduce  that  for  l^(z)  =  0, 
P{z,T)  =  P^(T)  and  for  E^(z)  =  1,  P{z,T)  =  P^T).  As  a 
result,  l^(z)  e  [0,1]. 

Two  Scenarios  of  the  grading  of  the  two  basic  component 
phases,  ceramic  and  metal,  through  the  wall  thickness  are 
considered. 

Case  (1):  The  phases  vary  symmetrically  through  the  wall 
thickness,  in  the  sense  of  having  full  ceramic  at  the  outer 
surfaces  of  the  plate  and  tending  toward  full  metal  at  the 
mid-surface.  Eor  this  case,  V^(z)  can  be  expressed  as 


N  is  termed  the  volume  fraction  index  which  provides  the 
material  variation  profile  through  the  plate  wall  thickness, 
(0  <  Af  <  00).  A  pictorial  representation  of  the  distribution 
of  the  constituent  materials  are  shown  below  in  figure  2. 


z 


Figure  2:  Distribution  of  the  constituent  materials 
through  the  plate  thickness  for  the  symmetric  case. 


Case  (2):  The  phases  vary  non-symmetrically  through  the 
wall  thickness,  and  in  this  case  there  is  full  ceramic  at  the 
outer  surface  of  the  plate  wall  and  full  metal  at  its  inner 
surface.  Eor  this  case,  V[,(z)  can  be  expressed  as 


Below  is  a  pictorial  representation  of  the  antisymmetric  case 
shown  in  figure  3. 


z 


Figure  3:  Distribution  of  the  constituent  materials 
through  the  plate  thickness  for  the  antisymmetric  case. 


yA^) 


l  +  sgn(z)^  ^  f l-sgn(z)^ 

h/2j  [  2  J  [h/2)  [  2  j 


(5) 


Where  the  signum  function  is  defined  as 


1,  z  >  0 

sgn(z)  =  <  0,  z  =  0  (6) 

-1,  z<0 


It  should  be  noted  that  in  contrast  to  case  (2),  where  there 
exists  coupling  between  stretching  and  bending,  such 
coupling  is  not  present  for  the  symmetric  case  (1).  Also,  for 
the  purposes  of  simplicity  the  Poisson’s  ratio  will  be 
assumed  to  be  constant  throughout  the  plate  structure.  From 
equations  (1)  -  (7),  the  effective  material  properties  of  a  EG 
plate  can  be  expressed  for  the  antisymmetric  case  as 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 

Page  2  of  1 1 


Proceedings  of  the  2009  Ground  Vehicle  Systems  Engineering  and  Technology  Symposium  (GVSETS) 


[e(z,  t),  p{z,  t)]  = 

{Em(T\Pm(T)) 


h  +  2z 
2h 


N 


(8) 


And  for  the  symmetric  case  as 


[E{z,T),p{z,T)]={E,^{np,^{T)) 


l  +  sgn(z 


/  A 

-z 

^f  l-sgn(z)^ 

[h/2j 

1  2  J 

+  {e^(T),  p^(T)) 

(9) 

v(z,r)  =  v(7’)  (10) 

Where, 

Ecm  ~  Ef-  —  Eyyi  ,  ~  Pc  ~  Pm  (  ^  l  ) 

3.  KINEMATIC  EQUATIONS 

3.1  The  3-D  Displacement  Field 

Consistent  with  the  classical  plate  theory  [6],  the 
distribution  of  the  3-D  displacement  quantities  through  the 
wall  thickness  can  be  expressed  as 

/  \  /  \  dwAx,  y,  t) 

u\x,  y,  z,t)  =  UQ\x,  y,  t)-z - , 

dx 

(  \  (  \  dwn,{x,y,t) 

v{x,y,z,t)  =  VQ{x,y,t)-z - - - ,  (12a-c) 

dy 

w{x,y,  z,t)=WQ{x,y,t) 


Within  these  equations,  (u,  v,  w)  are  the  3-D  displacement 
quantities  along  the  (x,  y,  z)  directions,  respectively.  While, 
Uo,r’o,andwo  are  the  2-D  displacement  quantities  of  the 
points  on  the  mid-surface. 

3.2  Non-Linear  Strain-Displacement  Relationships 

The  nonlinear  strain  displacement  relationships  across  the 
plate  thickness  at  a  distance  from  the  mid-surface  take  the 
form  [2,  3,  6] 


du  if  dx  1 

dx  2  fdx  J  dy  2 


dy 


-  du  dv  dw  dw 

2^  —  Y  — - 1 - 1 - ,  6  y  —  6 

xy  /  xy  ^  ^  ^  ^  ’  A4 

dy  dx  dx  dy 


yz 


=  0 


(13) 


Substitution  of  Equations  (12a-c)  into  Equations  (13)  results 
in  the  strain  measures  across  the  plate  thickness  in  terms  of 
the  2-D  displacement  quantities  of  the  mid-surface  of  the 
plate  expressed  as 


(14) 
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In  the  above  expressions,  £yy, K^y),  are  referred  to  as 
the  membrane  strains  and  fsxx >  ^yy >  Yxy')  referred  to  as 
the  flexural  bending  strains  which  are  also  known  as  the 
curvatures. 

4.  CONSTITUTIVE  EQUATIONS 

The  stress-strain  relationships  for  a  state  of  plane  stress  is 
expressed  as  [10] 


^xx 

"Gii 

Gi2 

0  ■ 
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The  material  stiffnesses,  Qij  (z),  (i  =  1,2,6)  are  given  by  [9, 

10] 


E, 


=  F  /)' 


1 

N  +  3 


— 1 
N +  2  4N +  4  J 


F 

12 


(21) 


2i  1  -  Q22  - 


e{z,t)  ^  vEizJ) 
- Y  ’  2i2  - - — 


(16) 


And  for  the  symmetric  case, 


Ei  = 


F  h 
N  +  l 


+  Eh,  E,  =0,  £3 


4(v  +  3)  12 


(22) 


The  standard  force  and  moment  resultants  of  a  plate  are 
defined  as 


,  ij  iJ 


( i,  j  =  x,y,xy) 


(17) 


With  the  use  of  Equations  (13)  -(17),  the  stress  resultants 
and  stress  couples  are  related  to  the  strains  by  [3] 


~[A]  [5]1p°^| 

[B]  [£>%(!)  I 


(18) 


In  which  [A],  [5],  and  [D]  are  the  respective  in-surface, 
bending-stretching  coupling,  and  bending  stiffnesses.  Eor  the 
case  of  symmetric  EG  Plates,  [E]  =  0,  since  there  is  no 
bending-stretching  coupling.  The  global  stiffness  quantities, 
Aij ,  Bij ,  and  Dij ,  (/,  )=1 ,2,6)  are  defined  as 


5.  GOVERNING  EQUATIONS 

Hamilton’s  principle  is  used  to  derive  the  equations  of 
motion  and  the  boundary  conditions.  It  is  formulated  as 

dJ  =5\{u  +V-K)dt  =  0,  (23) 

Jfo 


Where  tg  and  t^are  two  arbitrary  instants  of  time.  IJ  denotes 
the  strain  energy,  V  denotes  the  work  done  by  surface 
tractions,  edge  loads,  and  body  forces,  and  K  denotes  the 
kinetic  energy  of  the  3-D  body  of  the  structure,  while  8  is 
the  variational  operator.  In  Equation  (23),  the  strain  energy 
is  given  by 


5U  = 


1. 1 


+c>yySeyy+c7^dy^)dz 


dCln 


(24) 


where  fig  denotes  the  mid-surface  area  of  the  panel.  The 
work  done  by  external  loads  is  expressed  as 


In  view  of  Equations  (8),  (9),  (10),  (16),  and  (19),  the 
global  stiffness  quantities  are  determined  as 

[("^1 1  ’  "^22 )’  (^1 1  ’  ^22  ),  {Di  1 ,  D22 )]  =  -  Y  {Ei  ,  £2  >  ^3 )  2 

l-v 

{'^l2’Ei2,Di2)=  ^(£j,£2,£3), 

l-v 

("^66  ’  ^66  ’  ^66  )  =  Tf.  \  (^1  ’  ^2  ’  ^3  )  (20a-c) 

2(1 +  vj 


where  for  the  antisymmetric  case, 
E,.^h 


£1  = 


N  +  l 


_  E  E 


sv  =  -\ 


j*  Pt  {x,  y)Sw^x,  y,  ^  -  Ph  {x,  y)Sw^x,  y-  ^ j 


d^r 


CT^^du  -I-  cr^y(5v 


(25) 


In  the  above  expression,  Pf-{x,y')  is  the  distributed  force  at 
the  top  surface  (z  =  h/2).  Pi,  (x,y)  is  the  distributed  force  at 
the  bottom  surface  (z  =  —h/2),  ,  ayy  a*y ,  are  the 

specified  stress  components  along  the  plate  edges,  and 
(5u,  Sv)  are  the  virtual  displacements  along  the  normal  and 
tangential  directions,  respectively,  along  the  plate  edges. 
Considering  only  the  transversal  inertia  of  the  structure,  the 
kinetic  energy  is  given  by 
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(26) 


Which  implies  that  the  the  variation  in  kinetic  energy  is 
expressed  as 


dM  dM  gvvn  dwn 

V  = _ —  +  2 _ —  +  N  -  +  N  ^ 


dy 


dy 


(29b) 


h 

^0=1^^  >0(z)* 


I 


(30) 


h 

SK  =  ^  p{z)VSWdzdna  (27) 

where  p(z)  is  the  mass  per  unit  volume. 

Considering  Equation  (23)  in  conjunction  with  Equations 
(24)-(27)  ,  along  with  the  constitutive  equations  (15),  the 
strain-displacement  relationships,  equations  (14),  and 
carrying  out  the  integration  through  the  thickness, 
integrating  by  parts  whenever  feasible,  using  the  expression 
of  global  stress  resultants  and  stress  couples,  while  retaining 
only  the  transversal  load,  transverse  inertia,  and  transverse 
damping  results  in 


Invoking  the  arbitrary  and  independent  character  of 
variations  5uo,  ^Wg,  5(dwo/9x)  ,  and  5(9wo/dy)one 
obtains  the  equations  of  motion  and  as  a  by-product  the 
boundary  terms  or  conditions.  This  results  in  three  equations 
of  motion  in  terms  of  the  global  stress  resultants  and  stress 
couples  and  four  boundary  conditions  along  the  plate  edges. 
These  equations  of  motion  and  boundary  conditions  can  be 
respectively  expressed  as 
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(28) 


In  the  above  equation,  (t7j.  t2y),  and /g  (referred  to  as  the 
inertia)  are  given  by 


dM,,  dMy 
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dx 
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The  boundary  conditions  become 
Along  the  edges  x  —  0,Li 
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or 

o 

II 

(34b) 
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XX  XX 
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or 
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creates  a  spherical  shock  wave  which  travels  radially 
outward  in  all  directions  with  diminishing  velocity.  The 
form  of  the  incident  blast  wave  from  a  spherical  charge  is 
shown  in  figure  4.  Where  Pso  is  the  peak  overpressure  above 
ambient  pressure,  Pq  is  the  ambient  pressure,  is  the  time 
of  arrival,  tp  is  the  positive  phase  duration  of  the  blast  wave, 
and  t  is  the  time.  The  waveform  shown  in  figure  4  is 


Figure  4:  Incident  pressure  profile  of  a  blast 
wave. 


given  by  an  expression  Known  as  the  Eriedlander  equation 
and  is  given  as 


Eor  the  case  of  all  edges  simply  supported  and  freely 
movable  the  boundary  conditions  are  as  follows: 

^0  =N^  =0,N^  =Nl,  on  x  =  0,L^ 

Wo=M^=Ny,=0,N^=N‘;y  on  y^0,L^  (40) 

It  should  be  mentioned  for  clarification  sake  that  for 
compressive  edge  loading  N*^.  —  —N^^,  and  N*y  —  —Nyy . 

6.  AIR-BLAST  LOADING 


\  -« 


_ c 

t 

V  P  J 


Where 


Pso  - 


1772  _  114  108 

73  72  ^  Z 


(41) 


(42) 


With  the  ever  increasing  demands  for  increased  safety  for 
the  solders  in  the  field  to  operate  structurally  sound  vehicles 
in  the  event  of  an  lED  or  some  other  type  of  explosive,  it  is 
imperative  that  an  understanding  of  the  structural  response 
of  various  components  within  military  combat  vehicles 
under  an  explosive  blast  be  understood  so  that  measures  can 
be  taken  from  a  design  standpoint  to  ensure  the  durability 
and  survivability  of  these  components.  To  begin  to  achieve 
this  understanding,  the  type  of  explosive  loading  considered 
here  is  a  free  in-air  spherical  air  burst.  Such  an  explosion 


In  equations  (41)  and  (42)  Z  is  known  as  the  scaled 
distance  given  by  Z  =  with  R  being  the  standoff 

distance  in  meters  and  W  being  the  equivalent  charge 
weight  of  TNT  in  terms  of  kilograms.  Also,  a  is  known  as 
the  decay  parameter  which  is  determined  by  adjustment  to  a 
pressure  curve  from  a  blast  test. 

Eor  the  conditions  of  standard  temperature  and  pressure 
(STP)  at  sea  level,  the  time  of  arrival  t^,  and  the  positive 
phase  duration  tp ,  can  be  determined  from  [4] 
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t  _  R 

h  R\  v^l  y 


(43) 


Where  represents  either  the  arrival  time  or  positive  phase 
duration  for  a  reference  explosion  of  charge  weight  W^,  and 
t  represents  either  the  arrival  time  or  positive  phase  duration 
for  any  explosion  of  charge  weight  W.  The  determination  of 
the  standoff  distance  for  any  charge  weight  W  follows  a 
similar  reasoning.  The  application  of  these  relationships  is 
known  as  cube  root  scaling.  It  should  be  understood  that  in 
applying  these  relationships  that  the  standoff  distances  are 
themselves  scaled  according  to  the  cube  root  law. 

7.  SOLUTION  METHODOLOGY 

To  satisfy  the  first  two  equations  of  motion,  equations 
(31a,b),  a  stress  potential  will  be  utilized  which  allows  the 
in-plane  stress  resultants  to  be  expressed  by  letting 


D  = 


El  £3  E2 

E,(1-v2)  ’ 


El 


(47b, c) 


This  gives  one  governing  equation  with  two  unknowns, 
Wg  and  0.  One  more  equation  is  needed  in  terms  of  the  same 
unknowns  which  will  give  two  equations  in  terms  of  two 
unknowns  which  can  then  be  solved.  This  will  come  from 
the  compatibility  equation.  By  eliminating  the  in-plane 
displacements  from  the  strain-displacement  relationships, 
equations  (14)  the  relationship  between  the  in-plane  strains 
and  the  transversal  deflection  known  as  the  compatibility 
equation  can  be  shown  to  be  given  by 


^yy 

dx^ 


dY 


dxdy 


^  _ 


7 -,2 

O  Wq 

,  dxdy ^ 


dx^  d^y 


(48) 


In  view  of  equations  (44)  ,  (45),  and  (46),  the  compatibility 
equation  is  expressed  as 


N 


XX 


N 

dx^’  ^  dxdy' 


(44) 


The  third  equation  of  motion,  equation  (31c),  can  be 
expressed  in  terms  of  two  unknown  variables,  the  stress 
potential  0  and  the  transverse  displacement  Wg.  To 
accomplish  this,  a  partial  inversion,  of  equation  (18),  the 
details  of  which  are  not  presented  here,  needs  to  be  carried 
out.  Performing  a  partial  inversion  results  in  [4] 


"  [A*] 

rlBY  1  ^,3^ 


Where, 

kl-pp.  kl=4Ar'[si  .-kr-MAp 

[D*]=[Dl-[illA]A[il]  (46) 


Using  Equations  (44),  (45),  and  (46)  and  simplifying. 
Equation  (31c)  takes  the  form 


V  V  =  El 


/  9  \ 

D  Wq 

^  dxdy 


2  2 
d  Wq  d  Wq 

dx^  d^y 


-(l-v)V^V^ 


(49) 


In  equations  (47)  and  (49),  V^=  d^ldx^  +d^ldy^  where 
V  is  referred  to  as  the  Laplacian  operator. 

Equations  (47a)  and  (49)  are  the  basic  governing 
equations  used  to  investigate  the  structural  response  of  EG 
plates  under  external  excitation  loading.  Eor  the  purposes  of 
this  paper,  from  this  point  forward,  the  thermal  terms  will  be 
discarded.  To  this  end,  to  solve  equations  (47a)  and  (49),  the 
approach  adopted  from  [2]  will  be  utilized.  In  this  respect, 
the  following  functional  forms  are  assumed  for  Wg  and  0 
[2]. 


Wq  (x,  y ,  f )  =  (f )  sin  sin  y 


(50a,b) 


4>  =  A^n  (f)cos  2/l„x  -I-  (f  )cos  2//„y  -1- 

^mn  (0  cos2yl„xcos2//„y-h£)„„  (0sin2A„  xsin2//„j  + 


1  -KT*  2  ,  1  -KT*  2 

-N^y  +-NyyX 


DVY  - 


'dY_dY 

dy^ 


^  d'^cj)  ^ 
dxdy  dxdy 


5^0  d^WQ 
dy^  dx^  ^ 


IoWo+Cwo=Pt+0^^N'^ 


(47a) 


Where, 


Where  —  nn/L2^  m,n  —  1,2, ...  are  the 

number  of  half  waves  in  the  x  and  y  directions,  respectively, 
and  it)  is  the  amplitude  of  deflection.  Also, 

Aynn  (0<  ^mn  Cmn  (0<  (t)  are  coefficients  to  be 

determined.  By  substituting  equations  (50a,b)  into  equation 
(49),  the  coefficients  it),  Bmn  (t),  Cmn  (t),  and  it) 
are  determined  as 
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^mn  [t)  = - — G - [t)= — — 

32^2  32mI 

C„.„it)=D^^{t)  =  0 


(51) 


By  letting, 


Pf  (f)  =  Pf„„  (t)sin  sin  A^^xsin  y  (52a) 


xp  —  L1/L2  is  referred  to  as  the  aspect  ratio  and  O  = 
Nyy/Nxx  is  referred  to  as  the  compressive/tensile  edge  load 
ratio. 

Equation  (53)  is  a  nonlinear  equation  in  terms  of  the  plate 
deflections  as  a  function  of  time.  It  is  interesting  to  note  that 
equation  (53)  is  very  similiar  to  Duffing’s  Equation.  To 
obtain  the  plate  deflections  as  a  function  of  time,  equation 
(53)  is  solved  using  the  Eourth-Order  Runge-Kutta  Method. 


which  implies  through  integration  of  both  sides  over  the 
plate  area  that 


4  fG  rh 

Pmn  (0  =  y-j-  \PAt)  sin  sin  ydxdy  (52b) 

E\E2  Jo  Jo 


8.  RESULTS  AND  DISCUSSION 

To  validate  the  present  theory,  comparisons  are  made  with 
Akay  [1]  who  considered  a  step  loading  excitation  of  a 
simply  supported  elastic  plate  based  on  a  von-karman 
nonlinear  theory.  To  make  this  comparison  the  external 
excitation,  (t)  ,  expressed  as 


or  through  integration  gives 


P,(0  =  P[//(f)-//(f-fo)] 


Pmn  (0  = 

K 

and  introduction  of  equations  (50a,b)  and  (51)  into  equation 
(47a)  and  retaining  the  resulting  equation  along  with  the 
unsatisfied  boundary  conditions  in  the  energy  functional  and 
applying  the  Extended  Galerkin  technique  results  in  the 
following  nonlinear  differential  equation  governing  the 
structural  response  of  EG  plates,  under  external  excitation. 

{tWmnit)  ^33^ 

=  Pmn(t) 


Where,  w^n  (t)  represents  the  amplitude  of  deflection  of  the 
plate  as  a  function  of  time,  P^n  (0  =  16Pt(t)//o7r^,  oj^n  = 
V^mn  / 4  is  ths  natural  frequency  of  the  EG  plate,  and 

=  C /ZIqCd^j^  is  the  non-dimensional  damping  factor, 
and  +  Pn)/16  •  It  should  be  noted  that  at  the 

center  of  the  plate  (x,  y)  =  (Li/2  ,L2/2),  w^n  (t)  is  equal 
to  the  maximum  deflection  of  the  plate.  In  these  latter 
expressions, 

Kmn  =  +  2m^n +  « V* 


Where, 


was  applied.  P[{t)  is  referred  to  as  the  Heaviside  Step 
function  defined  as  //(t)  =  1  for  t  >  0  and//(t)  = 
0  for  t  <  0. 


Fig  5.  Comparisons  of  solutions  of  the  time  history  of  the  Central  deflection  Under 
a  Step  Load 


The  geometrical  and  material  properties  used  were, 

L]  =  2.438  m,  h  =  0.00635  m,  y/  =  l 
E^=103GPa,  p^=25Alkglm\  i^av.=0.25 
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In  addition,  the  constant  pressure  and  time  duration  were 
taken  as  E  =  48.82  Ea,  0.2  sec.  Since  only  a  metallic 

isotropic  plate  was  made  for  comparison,  the  volume 
fraction  index  was  taken  as  N  —  2000  which  implies  fully 
metal.  With  this  in  hand,  the  central  deflection  time  history 
is  displayed  in  Eig  5.  The  results  in  this  figure  reveal  close 
agreement  between  Akay’s  finite  element  method  and  the 
present  analytical/  approximation  method  employed  here  in 
the  present  analysis. 

To  illustrate  the  present  approach,  a  ceramic-metal 
functionally  graded  plate  consisting  of  Ti-6A1-4V  and 
Aluminum  Oxide  with  the  following  material  properties, 
which  were  adopted  from  [9],  were  considered  for  the 
numerical  results  presented. 

=  320.24  GPa,  /?,  =  3750  kg  /  =  0.26 

=105.7  GPa,  =4429kg/m^  =0.2981 

=0.2791 

The  geometrical  properties  used  for  the  EG  Plate  are 
L,  =lm,  ^//'  =  Lj/L2  =1,  and  unless  otherwise  stated 
/j  =  0.0254  m  and  the  halfwaves,  (m,n)  =  (1,1).  In  addition, 
the  following  reference  values,  in  Table  1.  were  utilized  to 
determine  the  time  of  arrival  and  positive  phase  duration  [5]. 


Table  1.  Airblast  Parameters  Versus  Distance 
for  a  One  Kilogram  (Wj )  TNT  Spherical  Air 
Burst  [5]. _ 


Standoff 

Distance, 

(m) 

Arrival 

Time, 

(msec) 

Positive  Phase 
Duration,  tp^ 
(msec) 

1.0 

0.532 

1.79 

In  Pig  6,  comparisons  of  the  central  deflection  of  the  plate 
for  various  voulme  fraction  indexes  is  depicted,  it  can  be 
seen  that  the  symmteric  functionally  graded  case  with  a 
volume  fraction  index  of  N=0.5  gives  the  lowest  deflections 
as  a  function  of  time,  as  compared  to  the  fully  metal 
isotropic  case. 

Pig  7.  Is  the  counter  part  of  Pig  6.  With  the  exception  of 
the  effect  of  damping.  The  effect  of  damping  shows  a  rapid 
attentuation  of  the  deflections  over  a  very  short  period  of 
time.  It  is  shown  that  damping  plays  a  major  role  in  the 
decrease  of  deflections  of  the  plate. 

In  Pig  8.  It  can  be  seen  that  the  effects  of  various  amounts 
of  damping  on  the  central  deflections  vs  time  attenuates 
faster  as  the  amount  of  damping  is  increased  for  a  fixed 
volume  fraction  of  the  constituent  materials.  The  effects  of 
the  compressive/tensile  edge  loading  on  a  fg  plate  for  a  fixed 
volume  fraction  index  and  a  fixed  amount  of  damping  on  the 


deflection-time  response  is  depicted  in  Pig.  9.  It  is  clear  that, 
for  tensile  edge  loading,  the  magnitudes  of  the  deflections 
decreases  in  contrast  to  compressive  tensile  edge  loading 
where  the  deflections  are  increased. 


Fig  6.  Implications  of  the  volume  fraction  index  on  the  central  deflection  as  a  function 
of  time.  {/^=0.766  m,  f^=0.000408s,  f^,=0. 001372s) 


Fig.  7.  The  counterpart  of  Fig  6.  with  a  damping  measure  of  A,  ,=0. 1 

Pig  10.  reveals  the  effects  of  the  aspect  ratio  on  the 
deflection-time  response  which  indicates  that  the  deflections 
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are  larger  for  smaller  aspect  ratios  than  for  larger  aspect 
ratios.  Erom  a  design  standpoint,  utilizing  larger  aspect  ratio 
panels  would  be  more  beneficial  by  give  smaller  deflections 
under  an  explosive-type  excitation,  hence  decreased  stresses. 


Kig  8.  Implications  of  various  amounts  of  damping  on  the  deflection-time  response 
of  an  antisymmetric  KG  Plate.  {N=0.5) 


Fig.  9.  The  implications  of  the  comprcssivc/tcnsilc  edge  loading  on  the  dcflcction-timc 
response  of  an  antisymmetric  FG  plate.  (A,  ,=0.2,  N=0.5) 


Fig.  10.  The  effects  of  various  aspect  ratios  of  an  antisymmetric  FG  plate  on  the 
deflection-time  response.  (N=0.5) 


9.  CONCLUDING  REMARKS 

A  rigorous  treatment  of  functionally  graded  plates  with 
grading  in  the  transverse  direction  has  been  studied. 
Validations  with  a  simpler  transversal  excitation  step  load 
found  in  Akay  [1],  who  took  a  specialized  finite  element 
approach,  has  been  made.  It  has  been  shown  that  damping 
has  an  important  effect  when  it  comes  to  the  attenuation  of 
the  deflections.  It  has  also  been  shown  that  other  factors 
such  as  the  compressive/  tensile  edge  loading,  the  aspect 
ratio,  and  the  symmetry  of  the  transverse  grading  throughout 
the  structure  plays  an  important  role  in  the  deflection-time 
history  of  the  structure. 

In  Eig.  6.  It  was  shown  that  functionally  grading 
inherently  reduces  the  deflections  when  compared  to  the 
isotropic  metallic  case.  Erom  a  design  standpoint,  it  would 
be  appropriate  at  this  point  to  state  that  integration  of 
functionally  graded  materials  within  plate-type  structures 
would  benefit  the  structural  response  of  the  structure.  Also, 
it  should  be  mentioned,  although  not  shown  here,  that  the 
choice  of  the  ceramic  and  metal  constituent  materials  chosen 
would  also  have  a  great  impact  on  the  response  of  the 
structure. 

The  idea  is  to  reduce  the  stresses  within  the  structure 
concerned  here.  By  reducing  the  magnitude  of  the 
deflections,  the  stresses  are  reduced.  It  is  hoped  and  realized 
that  this  present  study  presented  here  will  give  insight  into 
some  of  the  factors  that  can  play  an  important  role  in  the 
structural  response  of  functionally  graded  plates  and  fill  in 
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some  of  the  fundamental  missing  gaps  within  this  subject 

area. 
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